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1. Introduction 

When cosmic rays collide with the nuclei of the atmosphere, they initiate extensive air show¬ 
ers (EAS). The interpretation of EAS data relies to a large extent on the understanding of these 
air showers, specifically on the correct modeling of hadron-air interactions that occur during the 
shower development. Experiments such as the Pierre Auger Observatory [9], KASCADE-Grande [10], 
IceTop [11] or the Telescope Array [12] use models for the interpretation of measurements. How¬ 
ever, there is mounting evidence that current models give a poor description of muon production in 
air showers (see Refs. [13, 14, 15]). 

Unfortunately, there exist no comprehensive and precise particle production measurements 
for the most numerous projectile in air showers, the Ti-meson. Therefore, new data with pion 
beams at 158 and 350GeV/c on a thin carbon target (as a proxy for nitrogen) were collected by the 
NA61/SHINE experiment at the CERN SPS. 

Preliminary spectra of unidentified hadrons have been previously derived from fhis dafa sef and 
Ihe specfra revealed discrepancies befween fhe dafa and predictions from generafors for hadronic 
inferaclions [16, 17, 18]. 

In fhis confribufion we will presenf fhe measuremenf of specfra of idenfified charged pions and 
p° mesons in n^+C inferacfions af 158 and 350GeV/c. 

2. The NA61/SHINE Experiment 

NA61/SHINE* [1] is a multi-purpose fixed fargef experimenf fo sfudy hadron production in 
hadron-nucleus and nucleus-nucleus collisions af fhe CERN Super Proton Synchrofron (SPS). 
Among ifs physics goals are precise hadron producfion measuremenfs for improving calculafions 
of fhe neufrino beam flux in fhe T2K and Eermilab neufrino oscillation experimenf [2] as well as 
for more reliable simulafions of hadronic inferacfions in air showers. Moreover, p-i-p, p-i-Pb and 
nucleus-i-nucleus collisions are measured fo sfudy fhe properfies of fhe onsef of de-confinemenf and 
search for fhe crifical poinf of sfrongly inferacfing maffer (see e.g. Ref. [3]). 

The NA61/SHINE Collaboration uses large lime-projection-chambers (TPCs) inherifed from 
fhe NA49 experimenf [4] fo measure fhe charge and momenlum of parlicles. The momenlum 
resolution, a{l/p) = o{p)/p^, is abouf 10^^(GeV/c)^* af full magnetic field and fhe fracking 
efficiency is better lhan 95%. A sef of scinfillalion and Cherenkov counfers as well as beam po- 
sifion defecfors upsfream of fhe specfromefer provide timing reference, idenfificafion and position 
measuremenfs of fhe incoming beam parficles. Parlicle idenfificafion is achieved by measuring fhe 
energy loss along fhe fracks in fhe TPCs and by defermining Iheir velocify from fhe lime of llighl 
provided by large scinlillalor walls placed downslream of fhe TPCs. The cenlralily of nucleus- 
nucleus collisions can be estimated using fhe measuremenf of fhe energy of projecfile speclafors 
wilh a calorimeler [5] localed behind fhe time of flighl defecfors. Eor nucleon-nucleus collisions, 
fhe cenlralily is determined by counting low momenlum parlicles from fhe large! (so called ‘gray 
protons’) wilh a small TPC around fhe large! [6]. 

Dafa faking wilh fhe NA61/SHINE experimenf slarled in 2007. Afler a firsl run wilh prolon 
on carbon af 31 GeV/c, fhe dafa acquisifion syslem was upgraded during 2008 fo increase fhe evenl 

'shine = SPS Heavy Ion and Neutrino Experiment 
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recording rate by a factor of sslO [7]. In the last years, a wealth of data has been recorded by 
the experiment at beam momenta ranging from 13 to 350GeV/c with various beam particles and 
targets. In this paper we present results obtained from a special run with negative pions as beam 
particles and carbon as the target. Since pions are the most numerous particles in an air shower, 
this data will help to improve the interpretation of air showers at ultra-high energies. 


3. Spectra of Charged Pions 

For each track detected in the TPCs of NA61, the particle type can be estimated by using 
the truncated mean of the energy that is deposited per unit track length (dE/Ax) along the particle 
trajectory. An example of a dF/dv-distribution in a specific bin in momentum p and transverse 
momentum pj is shown in Fig. 1. As can be seen, the distribution can be well described by the 
sum of the energy loss distributions of electrons, protons, pions and kaons (see Ref. [8] for details) 
and given the fitted fraction of each particle type, the corresponding number. An, of produced tracks 
within each p /px-bin can be reconstructed. 

This number is then corrected for 
the detector acceptance, selection ef¬ 
ficiency, feed-down from weak de¬ 
cays and re-inferacfions in fhe far- 
gef. The laffer fwo correcfions are cur- 
renfly esfimafed using model predic¬ 
tions (Epos 1.99, QGSJetII-04, DP- 
MJet 3.06) and fhey are fypically well 
below 5%, buf can reach up fo 20% 
af low particle momenfa. Overall, fhe 
sysfemafic uncertainly of fhe correcled 
number of hacks. An', is esfimafed fo 
be <7%. 

The average mulfiplicify of parli- 
cles produced wilhin a p/pj-bin is fhen 
obfained by dividing An' by fhe folal 
number of evenfs in which an inferac- 
fion occurred, Aprod- Aprod is esfimafed 
by exfrapolaling fhe number of recorded 
inleraclion higgers fo full phase space. 

The correction amounfs fo (92.5 ± 3.5)% al 158 GeV/c and (92.5 ±4.0)% al 350GeV/c, where fhe 
uncerlainfy was derived by running differenl generators fo evaluale fhe correction. 

The measured average mulliplicilies of charged pions are shown in Fig. 2 and fhe measure- 
menls are compared fo predictions of pion produclion in n^+C al 158 GeV/c from hadronic inler- 
acfion models in Fig. 3. As can be seen, none of fhe generalors describes fhe dala well. 







Figure 1: Example of a dE/dx fit ((p) = 28.4GeV/c, 0 < 
Pt < 0.1 GeV/c). The middle panel shows the energy de¬ 
posit of positively and negatively charged tracks. The fitted 
particles are indicated by colored histograms. Residuals to 
the fit are shown in the upper and lower panels. 
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Figure 2: Inclusive spectra of charged pions in 7r~+C interactions at beam energies of 158 and 350GeV/c. 
For better visibility, the spectra from the n'-^ momentum bin are multiplied by a factor of 1 / 4 ". The momentum 
increases from top to bottom as indicated in the legend on the right. Systematic errors are less than 7% 


4 . Spectra of p® Mesons 

The measurement of resonances in n+C is useful to constrain the production of meson, 
which is important to predict the number of muons observed in air showers as well as the baryon 
fraction (see e.g. Ref. [19]). 

In the inclusive invariant mass spectra there is a large combinatorial background, which 
dominates over the effective mass distributions of individual resonances. The method used to 
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Figure 3: Comparison of the measured production spectra of charged pions to predictions from hadronic 
interaction models used for the interpretation of cosmic ray data [25, 23, 24, 22]. The colors denote the ratio 
of data to simulation and the color scale is truncated at 0.5 and 2. The two empty p-bins at p < 2 GeV /c and 
p > 3 GeV/c are momenta excluded from the analysis due to the ambiguity in identification of pions. 
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Figure 4: invariant mass distribution in n^+C interactions at 158 GeV/c in the range 0.4 < .rp < 0.5. 

Dots with error bars denote the data and the htted resonance templates are shown as hlled histograms. The 
vertical lines indicate the range of the ht. 

estimate the background is the so called charge mixing, which uses the + 71^71^) mass 

spectra as an estimate of the background. 

The fitting procedure uses templates of the tt+tt^ mass distribution for each resonance. These 
templates are constructed by passing simulated tt+C interactions, generated with the EPOS 1.99 [20] 
hadronic interaction model using Crmc [21] (vl.5.3), through the full NA61 detector Monte Carlo 
chain. All the cuts that are applied to the data are also applied to the templates. This method of 
using templates allows for the fitting of both resonances with dominant three body decays, such 
as the ft), and resonances with non-TT+Ti^ decays, such as the K*^. The data is split into bins of 
Feynman-x, xp- 

The fit to the tt+tt^ mass spectrum is performed between masses of 0.4GeV/c and 1.5 GeV/c 
using the expression 


F{m) = 


where jS, is the relative contribution for each template, 7] , used. An example of one of these fits can 
be seen in Fig. 4, The templates in the fit are the background found from charge mixing and the 
following resonances: p°, K*°, (O, f 2 , fo (980), a 2 , T] and K°. 

The fitting method is validated by applying the same procedure to the simulated data set which 
was used to construct the templates for the fit. For the majority of xp bins there is good agreement 
between the fit and the true value, with some discrepancies for larger xp bins of up to 20%. This 
bias is corrected for in the final analysis. The dafa is also correcfed for losses due fo fhe accepfance 
of fhe defecfor, as well as any bias due fo fhe cufs used and any reconsfrucfion efficiencies. Aparf 
from fhe accepfance, fhese corrections are fypically less fhan 20%. 


The average mulfiplicify of p° mesons is presenfed in Fig. 5. Also shown are predicfions by 
EposI.99 [20], DPMJET3.06 [22], Sibyll2.1 [23], QGSJetII- 04 [24] and EposFHC [25]. 
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Figure 5: Average multiplicity of the meson in n+C at pbeam = 158GeV/c as a function of Feynman-x. 
The bars show the statistical errors; the bands indicate systematic errors. The lines depict predictions of 
hadronic interaction models: red - EPOS 1.99, blue - DPMJet 3.06, black - Sibyll 2. 1, green - QGSJetII- 
04, dashed red - EPOSLHC. 

It can be seen that there is an underestimation of the yield for almost all hadronic interaction 
models, with the exception of QGSJetII- 04 for xp > 0.8. It is interesting to note that while 
QGSJetII- 04 and EPOSLHC were tuned to NA22 7r++p data [26], there is an underestimation in 
n^+C. 

Systematic errors are estimated by comparing correction factors for different hadronic interac¬ 
tion models (EPOS and DPM Jet), comparing the correction for the bias using different background 
estimates and varying the cuts applied to the data. The systematic is dominated by the background 
estimates, up to 14%, where as the other errors are less than 4%. Other sources of uncertainty, such 
as using templates from a different model, are found to be much smaller. 

5. Conclusions and Outlook 

In this article, we summarized results from pion-carbon interactions measured with the multi¬ 
purpose experiment NA61/SHINE at the CERN SPS, which are of importance for the modeling of 
cosmic ray air showers. 

The comparisons to hadronic interaction models shown in this article suggest that these models 
require further tuning to reproduce the charged pion and p° meson spectra. 

It is planned to further refine both analyses presented here, including the measurement of 
inclusive spectra of charged kaons and protons as well as the study of the multiplicities of other 
resonances in addition to the p°. 
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